if !d.name eq 'PS' then begin
    device,xsize=18,ysize=12,yoffset=3
    !p.thick=1 & !x.thick=1 & !y.thick=1
end
!x.margin=[10,3]
!p.charsize=1.5

@eu_setup
!p.charsize=2.5 
nt=(size(u))[3]
nth=(size(th))[1]
print,nt
nl=64
year = 365.*24.*60.*60.
time = (nxaver*dt*indgen(nt))/year
good=where(time gt 2)
ng=(size(good))[1]

umax=max(abs(bx[*,*,good]))
umin=-bmax
lev=grange(umin,umax,nl)
rlev=34
uphi=reform(u[*,rlev,good])
umean=total(uphi,2)/double(nt)
tosc=fltarr(nth,ng)
for i=0,ng-1 do begin
 tosc[*,i]=uphi[*,i]-umean[*]
endfor
lev = 2.*max(abs(tosc))*indgen(nl)/float(nl-1) $
        -max(abs(tosc))
contour,rotate(tosc,3),time[good],th,/fi,nl=64,ytitle='!8B!D!7u!N!6',ystyle=1,xstyle=1

!p.multi=0
end
